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Figure 1: Phase diagram in the (x, T) plane [|12|]; here x represents the doping and T the temperature. 



1. Introduction 

Quantum electrodynamics in 2+1 dimensions (QED3) is an interesting testfield for understand- 
ing the mechanism of confinement in gauge theories and for the effective description of low- 
dimensional, correlated, electronic condensed matter systems, like spin systems or high-r c 
superconductors [ph. The compact formulation of QED3 is more suitable for studying the mech- 
anism of confinement, while both compact ^ and noncompact formulations arise in condensed 
matter systems. This work is devoted to clarify the relationship between these two formulations of 
QED3 on the lattice (for more details, see Ref. [[j]]). 

Compact QED3 without fermion degrees of freedom is always confining [^] ; a charge and anti- 
charge pair is confined by a linear potential, as an effect of the proliferation of magnetic monopoles. 
If matter fields are introduced, the interaction between monopoles could turn from l/x to — ln(x) 
at large distances x [|7|], so that the deconnned phase may become stable. However, it has been 
proposed that compact QED3 with massless fermions is always in the confined phase 

At finite T noncompact QED 3 is relevant in the analysis of the pseudo gap phase [ ^Q| ] of 
cuprates. In Fig. [j] we report the phase diagram of high-r ( cuprates. The small-x phase (x is the 
doping) is characterized by an insulating antiferromagnet (AF); by increasing x, this phase 
evolves into a spin density wave (SDW), that is a weak antiferromagnet. The pseudo gap phase 
is located between this phase and the J-wave superconducting (dSC) one. The effective theory 
of the pseudo gap phase JlQt ] turns out to be QED3 [||, |l2|, with spatial anisotropics in the 
covariant derivatives JrT|], and with Fermionic matter given by spin- 1/2 chargeless excitations of 
the superconducting state (spinons). These excitations are minimally coupled to a massless gauge 
field, which arises from the fluctuating topological defects in the superconducting phase. The SDW 



order parameter is the chiral condensate (yn//) [13]. If (y/y/) is different from zero, then the <f-wave 



superconducting phase is connected to the spin density wave one (see Fig. [l] case b); otherwise the 
two phases are separated at T = by the pseudo gap phase (see Fig. [T] case a). 

The determination of the critical number of flavors, Nf tC , such that for Nf < Nf iC the chiral 
condensate is nonzero, is a crucial point addressed in many works (see Ref. [JTJ] for a detailed list). 
Here, we shall not try to ascertain Nf jC . However, we will handle the problem of confinement in 
presence of massless fermions, by looking at the relation between monopole density and fermion 
mass. Moreover we shall focus on a comparative study of the two lattice formulations of QED3, the 
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compact and the noncompact ones. In particular, we revisit the analysis of Fiebig and Woloshyn of 
Refs. [[If], [l7|], where the dynamic equivalence between the two formulations of (isotropic) QED3 
is claimed to be valid in the finite lattice regime. 

2. Compact and noncompact formulations 

We adopt here the definition of the QED3 parity-conserving continuum Lagrangian density in 



Minkowski metric given in Ref. [18], 



& = - + yjiZ> M 7^ Xjfi - moViVi , (2-1) 

where (z = 1, .. . ,Nf) are 4-component spinors. Since QED3 is a super-renormalizable theory, 
dim[e] = +1/2, the coupling does not display any energy dependence. For the definition of the 7 
matrices, for the chiral and parity properties of the theory we refer to Ref. 

The lattice Euclidean action [jl9|, [l4|] using staggered fermion fields ~XjX> is given by 



N 

S = S G + ^^XMWn.nXiH , (2-2) 

i=l n,m 



where So is the gauge field action and M n m is the fermion matrix. The action ( |2.2| ) allows to sim- 
ulate N = 1,2 flavours of staggered fermions corresponding to Nf = 2,4 flavours of 4-component 
fermions \ff [|20|]. For the compact formulation So is the standard Wilson action written in terms of 
the plaquette variable U^ v (n). Instead, in the noncompact formulation one has 

Sg[oc] = j £ F MV (n)F MV (n) , F^ v (n) = {a v (n + fl) - CC v (n)} - {a^(n + v) - a^(n)} , (2.3) 

where (n) is the phase of the link variable U^(n) and j8 = 1/ (e 2 a). 

Monopoles are detected in the lattice using the method given by DeGrand and Toussaint [|]]] . 
In the noncompact formulation of QED3 monopoles are not classical solutions as in compact QED3, 
but they could give a contribution to the Feynman path integral owing to the periodic structure of 



the fermionic sector []22[]. 

As a signal for continuum physics, we look for plateau of dimensionless observables, such 
as p 2 (xx). There are two regimes: for j8 larger than a certain value, the theory is in the con- 
tinuum limit, otherwise the system is in a phase with finite lattice spacing, describing a lattice 



condensed-matter-like system. In Refs. [16, 17] it is shown there that, when (xx) is plotted ver- 
sus the monopole density p m , data points for QED 3 with Nf = and Nf = 2 in the compact and 
noncompact formulations in the lattice regime fall on the same curve to a good approximation (see 
Fig. § (left)). This led the authors of Refs. [16, 17] to conclude that the physics of the chiral sym- 



metry breaking is the same in the two theories. We want to verify, by the same method, if this 
conclusion can be extended to the continuum limit. 

3. Numerical results 

Our Monte Carlo simulation code was based on the hybrid updating algorithm, with a micro- 
canonical time step set to dt = 0.02. We simulated one flavour of staggered fermions corresponding 
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Figure 2: Left: Correlation between (xx) an d Pm for the compact (circles) and the noncompact (boxes) 
theories for Nf — 2 and 8 3 lattice according to Ref. [ jl7| ] (left) and to our results (right). 

to two flavours of 4-component fermions. Most simulations were performed on a 12 3 lattice, for 
bare quark mass ranging in the interval am = 0.01 -^0.05. We made refreshments of the gauge 
(pseudofermion) fields every 7 (13) steps of the molecular dynamics. In order to reduce autocor- 
relation effects, "measurements" were taken every 50 steps. Data were analyzed by the jackknife 
method combined with binning. 

As a first step, we have reproduced the results by Fiebig and Woloshyn which are shown in 
Fig. ^| (left). We find that also in our case data points from the two formulations nicely overlap (see 
Fig. (right)). 

Then, in Fig. || (left) we plot data for fi 2 (~XX) obtained in the compact formulation versus 
j5m. We restrict our attention to the subset of j8 values for which data points fall approximately on 
the same curve (this indicating the onset of the continuum limit) which in the present case means 
j8 = 1.9,2.0,2.1, corresponding to L//3 = 6.31,6.00,5.71. The ratio L//3 fixes the physical volume 
that is pratically constant in the considered range. A linear fit of these data points gives # 2 /d.o.f. 
~ 8.4 and the extrapolated value for j3m — > turns out to be f3 2 (x~X) = (1.54 ±0.25) x 10~ 3 . 
Restricting the sample to the data at jS = 2.1, the # 2 /d.o.f. lowers to ~ 1.3 and the extrapolated 
value becomes p 2 (xx) = (0.94 ±0.28) x 10~ 3 , thus showing that there is a strong instability in 
the determination of the chiral limit. If instead a quadratic fit is used for the points obtained with 
j8 = 1.9,2.0,2.1, we get P 2 (xx) = (0.91 ±0.45) x 10~ 3 with £ 2 /d.o.f. ~ 8.7. Owing to the large 
uncertainty, this determination turns out to be compatible with both the previous ones. 

In Fig. ^ (right) we plot data for f) 2 (x~x) obtained in the noncompact formulation versus /3m. 
Following the same strategy outlined before, we restrict our analysis to the data obtained with 
/3 = 0.7,0.75,0.8, which correspond to L/j8 = 17.14, 16, 15. If we consider a linear fit of these 
data and extrapolate to j8m 0, we get p 2 (x~X) = (0.45 ±0.03) x 10~ 3 with # 2 /d.o.f. ~ 17. 
Performing the fit only on the data obtained with j3 = 0.8, for which a linear fit gives the best 
# 2 /d.o.f. value ~ 16, we obtain the extrapolated value f3 2 {xx) = (0.66±0.07) x 10~ 3 . Therefore, 
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Figure 3: j8 (XX) versus pm in the compact (left) and in the noncompact (right) formulation. 
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Figure 4: (xx) versus p,„ in both the compact and the noncompact formulations on a 12 3 lattice. 



also in the noncompact formulation the chiral extrapolation resulting from a linear fit is largely 
unstable. A quadratic fit in this case gives instead a negative value for p 2 (xx)- 

The comparison of the extrapolated value for p 2 (~XX) m the two formulations is difficult owing 
to the instabilities of the fits and to the low reliability of the linear fits, as suggested by the large 
values of the # 2 /d.o.f. Taking an optimistic point of view, one could say that the extrapolated 
P 2 (XX) f° r j3 =2.1 in the compact formulation is compatible with the extrapolated value obtained 
in the noncompact formulation for j8 = 0.8. 

It is worth mentioning that our results in the noncompact formulation are consistent with 
known results: indeed, if we carry out a linear fit of the data for j8 = 0.6,0.7,0.8 and am=0.02, 
0.03, 0.04, 0.05 and extrapolate, we get P 2 (~XX) = (1-30 ±0.07) x 10~ 3 with an admittedly large 
X 2 /d.o.f. ~ 20, but very much in agreement with the value fi 2 (~XX) = (1-40 ±0.16) x 10~ 3 ob- 
tained in Ref. [19]. We stress that our results are plagued by strong finite volume effects, therefore 
our conclusions on the extrapolated values of fi 2 (~XX) we significant only in the compact versus 
noncompact comparison we are interested in. 

In Fig. |] we plot (xx) versus the monopole density p m . Differently from Fig. |2[ it is not 
evident with the present results that the two formulations are equivalent also in the continuum 
limit, although such an equivalence cannot yet be excluded. We arrive at the same conclusion 
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considering a 32 3 lattice. 

We studied also /3 3 p m versus /3m for the two formulations (for plots and more details, see 
Ref. [P). Our result show that data at different /3 values do not fall on a single curve, this suggesting 
that the continuum limit has not been reached for the monopole density for j8 = 2.2 in the compact 
formulation and j3 = 0.9 in the noncompact one. We found, however, that the monopole density 
is independent from the fermion mass. Since the mechanism of confinement in the theory with 
infinitely massive fermions, i.e. in the pure gauge theory, is based on monopoles and since the 
monopole density is not affected by the fermion mass, we may conjecture that this same mechanism 
holds also in the chiral limit. 
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